Monte Carlo method for evaluation of surface emission rate measurement uncertainty 


Yuan-Qiao Li,!:* Min Lin,!* Li-Jun Xu,!? Rui Luo,}? Yu-He Zhang,'? Qian-Xi Ni, t and Yun-Tao Liu!? 


‘China Institute of Atomic Energy, Beijing 102413, China 
°National Key Laboratory for Metrology and Calibration Techniques, Beijing 102413, China 
Department of Radiation Oncology, Hunan Cancer Hospital, Changsha 410031, China 
“The Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha 410031, China 


The aim of this study is to evaluate the uncertainty of 27a and 276 surface emission rates using the win- 
dowless multiwire proportional counter method. This study used the Monte Carlo method (MCM) to validate 
the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method. A dead-time mea- 
surement model for the two-source method was established based on the characteristics of a single-channel 
measurement system, and the voltage threshold correction factor measurement function was indirectly obtained 
by fitting the threshold correction curve. The uncertainty in the surface emission rate was calculated using the 
GUM method and the law of propagation of uncertainty. The MCM provided clear definitions for each input 
quantity and its uncertainty distribution, and the simulation training was realized with a complete and complex 
mathematical model. The results of the surface emission rate uncertainty evaluation for four radioactive plane 
sources using both methods showed an uncertainty consistency E, < 0.070 for the comparison of each source, 
and the uncertainty results of the GUM were all lower than those of the MCM. However, the MCM has a more 
objective evaluation process and can serve as a validation tool for GUM results. 
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I. INTRODUCTION 


Spills, leaks from unsealed sources, or breakage or breach 
of the integrity of sealed sources can result in radioactive 
contamination of material surfaces. This requires the effec- 
tive monitoring of surface radioactive contamination in ra- 
diation workplaces. Large-area œ- and -plane sources are 
widely used to provide the necessary traceability for mea- 
surements [1-5]. For absolute measurements of œ and 8 
surface emission rates, large-area windowless multiwire gas- 
flow proportional counter systems are commonly used. In this 
approach, the measurement function of the surface emission 
rate Æ is modeled as follows [6-9]: 
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where R is the particle count rate, B is the background count 
rate, T is the dead time of the counting system, and f is the 
correction factor. R and B can be obtained by direct measure- 
ment. However, for the system dead time 7 and correction 
factor f, there are differences in the measurement models ob- 
tained using different measurement instruments and methods. 

The conventional approach to the evaluation of surface 
emission rate uncertainty is the Guide to the Expression of 
Uncertainty in Measurement (GUM) method [10], which uses 
a nonlinear functional model to propagate uncertainty. This 
method determines the combined standard uncertainty of the 
output quantities using the law of uncertainty propagation. 
In the measurement function model of the surface emission 
rate E, the standard uncertainties of R and B are obtained 
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by statistically determining the experimental standard devi- 
ation of the measured values, a process known as Type A 
evaluation of uncertainty. For counting measurement systems 
based on single-channel analysis techniques, the system dead 
time 7 is typically obtained using the two-source measure- 
ment model, and the correction factor f is typically referred 
to as the voltage threshold correction factor fa. The value 
of fa is obtained by fitting the threshold extrapolation curve 
using the low-energy correction method. These two evalu- 
ation processes are usually classified as Type B uncertainty 
evaluations. This is due to the fact that when evaluating the 
uncertainty of surface emission rate using the GUM method, 
calibrators often make artificial assumptions about the proba- 
bility distribution of measurement results from complex non- 
linear equation models, usually assuming a rectangular distri- 
bution, and the effect of this method on the uncertainty eval- 
uation results is unknown. 

Owing to the nonlinear nature of the measurement model 
and the challenges associated with calculating the partial 
derivatives of the model in the surface emission rate uncer- 
tainty, the uncertainty can be measured using the Monte Carlo 
method (MCM) proposed in ISO/IEC GUIDE 98-3 [11]. As 
shown in Fig. 1, the MCM obtains the probability distribu- 
tion function (PDF) of the output quantity Y by randomly 
sampling the PDF of the input quantity X;. On this basis, the 
best estimates of output quantity, standard uncertainty, and 
coverage interval are obtained. 

The MCM has been used less frequently in uncertainty as- 
sessment studies of nuclear measurements [12, 13]. This is 
due to the fact that, although the GUM method has a gen- 
eralized law of propagation of uncertainty, a complete and 
specific measurement function model may differ due to the 
variations in experimental measurements for each measured 
quantity and the differences in function models for each com- 
ponent, resulting in a complex law of propagation of uncer- 
tainty with Taylor expansions. The staff may intentionally 
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Fig. 1. Description of evaluation uncertainty by the Monte Carlo 
method (MCM). The probability distribution of the input quantity 
X; involves the normal distribution gxi(&1), rectangular distribu- 
tion gx2(€2), triangular distribution gx3(€3), among others, and the 
probability distribution gy (7) of the output Y is unknown. 


simplify the measurement model through known evaluation 
results, avoid complex terms in the evaluation process, or treat 
them as constants to regain uncertainty with a smaller value, 
which may affect the accuracy of the traceability of quantity 
values [14]. In contrast, the MCM can directly perform calcu- 
lations for a variety of complex models to obtain the best esti- 
mate and standard uncertainty without considering the trade- 
offs between higher-order terms of the Taylor expansion in 
the law of the propagation of uncertainty. Therefore, the key 
to the MCM is to obtain a model of the function between the 
input and output quantities and its probability distribution for 
each measured quantity. 

To evaluate the 27a and 27 surface emission rate un- 
certainties, this study first describes the experimental setup, 
measurement method, and model for the system dead time 
and voltage threshold correction factor. The MCM is used to 
refine the surface emission rate measurement function model 
based on Eq. 1. Finally, the law of propagation of the uncer- 
tainty of the GUM is verified and compared with the func- 
tional model of the MCM. 


Il. METHODS 
1. Experimental set-up 


In this study, a large-area windowless multiwire gas-flow 
proportional counter [15] developed by the China Institute of 
Atomic Energy (CIAE) was used as the detector. The work- 
ing principle of this detector is based on the primary ion- 
ization of a and 8 particles in the sensitive volume of the 
proportional counter, which produces positive ions and elec- 
trons. These particles are then accelerated by an electric field 
to gain sufficient energy to trigger secondary ionization with 
gas molecules. The number of electron-ion pairs produced 
by this process is proportional to the number of electron-ion 


pairs produced by the primary ionization, resulting in a linear 
amplification of the primary ionization pulse signals from the 
a and (3 particles. 

For the single-channel counting measurement require- 
ments of this experiment, the electronic equipment included 
a high-voltage DC power supply, current-type pre-amplifier, 
pulse-shaping circuit, timing single-channel analyzer, monos- 
table control unit, and timer/counter. A diagram of the in- 
strument is shown in Fig. 2. P-10 gas (a mixture of 90% Ar 
and 10% CH4) was used as the operating gas at a flow rate 
of 100 mL/min. The working high voltage for the a-plane 
source was 2100 V, whereas that for the 3-plane source was 
2800 V. Measurements were performed on four plane sources 
C4 Am, 288Pu, Sr, and 2°4T1). According to oscilloscope 
observations, the pulse width of the pulse signal output from 
the main amplifier was generally less than 2 us. The output 
pulse width of the single-channel analyzer was approximately 
650 ns, whereas the pulse width of the monostable control 
unit had an adjustable range of 0.4—4 us. 


2. Dead time measurement 


The surface emission rate of plane sources generally does 
not exceed 5000 s~!, resulting in a dead time that produces 
leakage counts that may have a small effect on the count rate 
measurements [16-22]. However, because the windowless 
proportional counter is an absolute measurement device, it is 
necessary to measure the dead time of the system to correct 
for leakage counts. The dead time provided by the monos- 
table control unit within the measurement system is a non- 
extended dead time; however, the system dead time is also 
affected by various components. Therefore, it is necessary 
to measure the system dead time experimentally. ICRU Re- 
port 52 [23] recommends three methods for measuring the 
dead time. In this experiment, a two-source two-pulse method 
was used to measure the dead time of the system. In this 
method, pulse signals are generated from two different ra- 
dioactive sources and input into the detector separately to ob- 
tain the measurement count rates of the two sources, nı and 
ng. The two sources are then measured simultaneously using 
a detector to obtain the mixed measurement count rate n12. 
Ideally, the real count rate of three sets of pulse count rates 
ni, nh, and n4 should satisfy n} + n4 = nio, and according 
to the formula for the real count rate of a non-extended dead 
time system, n’ = n/(1 — Tn), the dead time formula for the 
two-source method is finally obtained: 
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3. Threshold correction factor measurements 


To avoid the influence of noise on count measurement, a 
discriminating voltage threshold must be set. Because the en- 
ergy of a particles is high and the discriminating threshold 
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Fig. 2. Actual size of a well-prepared line drawing. 


is set low, small changes in the a-particle count rate mea- 
surement can be ignored. However, the energy spectrum of 
6 sources is continuously distributed, and it is difficult to 
distinguish their low-energy signal from noise signals using 
the pulse amplitude. Consequently, setting the discriminat- 
ing voltage threshold may result in a loss of the -particle 
count rate, and a threshold correction for 6 sources must be 
performed [5, 24-26]. 

Some National Metrology Laboratories [27, 28] that use 
multichannel analyzers to construct measurement systems 
typically calibrate the site corresponding to the lower thresh- 
old by using 10% of the site of the °°Fe full-energy peak. This 
is primarily designed to correct for count-rate variations ow- 
ing to fluctuations in the gain of the working gas. Although 
a measurement system based on a single-channel analyzer 
cannot set a lower threshold by determining the absolute en- 
ergy, it is possible to set different lower discriminating volt- 
age thresholds to obtain the count rate at different discrimi- 
nating voltage thresholds, and then extrapolate to obtain the 
count rate when the threshold is zero. 

The threshold extrapolation curves are influenced by low- 
energy noise interference and the characteristics of the beta 
smooth continuum spectrum and can only be assumed to be 
a straight line over a small threshold interval. As a rule of 
thumb, the normalized value (predicted value/fitted intercept) 
of the fitted curve should not exceed 90%. The threshold ex- 
trapolation curve fitted using the least-squares method [29- 
33] is a straight line, where x represents the voltage threshold 
and y represents the measured count rate. The equation for 
the fitted curve is as follows: 


y =pr+q, (3) 


where p and q represent the slope and intercept of the line, 
respectively. This equation is solved as follows: 
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where ¿ = 1,2,...,9. To fit a straight line, a minimum of 
j = 5 measurements must be substituted into the equation to 
complete the fit with degrees of freedom v = j — 2. 


4. Evaluation by the GUM method 


The evaluation process of the GUM method is divided into 
two types of evaluations: A and B. This study focuses on 
two uncertainty components: system dead time and voltage 
threshold correction factor, and the process is a Type B eval- 
uation. The uncertainty obtained from the Type B evaluation 
of uncertainty can be obtained using Eq. (6). 


up = a/k = max(Aa)/k, (6) 


where a represents the half-width of the interval of the 
possible measurement values, or the maximum fluctuation 
max(Aa) if the coverage interval distribution is asymmet- 
ric [6, 34], and & is the coverage factor. Therefore, obtaining 
the coverage interval and probability distribution of the com- 
ponents is crucial. 

a. Evaluation of system dead-time uncertainty In the 
GUM method, the monostable control unit produces a square 
wave as the output signal. The pulse width of this signal 
is typically measured using an oscilloscope, which requires 
debugging the oscilloscope’s resolving power and sampling 
time interval. After debugging, the system dead-time inter- 
val half-width a is determined based on the time interval of 
the index value, as shown on the oscilloscope. Because the 
probability distribution function of the system’s dead time is 
unknown, it is generally assumed that the error follows a rect- 
angular distribution, in which case k = V3. 

b. Evaluation of voltage threshold correction factors un- 
certainty In this study, the threshold correction curves were 
fitted using the least-squares method. The expected correc- 
tions obtained from the calibration curves and their standard 
uncertainties were calculated using their intercepts, slopes, 
estimated variances, and covariances. Multiple sets of thresh- 
olds and the corresponding measurements were selected and 


substituted into Eqs. 4 and 5 to obtain the p and q values, 
which were substituted into the following equation: 
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where s(p) is the standard deviation of the slope, s(q) is the 
standard deviation of the intercept, s? is the estimated vari- 
ance of the fitted curve, and r(p, q) is the correlation coeffi- 
cient between the p and q estimates. This provided the com- 
bined standard uncertainty [10] of the measured count rate (y) 
corresponding to the specified prediction threshold (x) on the 
fitted curve. 


us(y) = u?(q) +a *u?(p) + 2a *u(p) *u(q) *r(p,q) AD 


The coverage interval [35, 36] for each measured count rate 
on the fitted curve was obtained as [px; + q + Uc(yi), puis + 


q + ue(yi)]. In addition, because y; is an experimental mea- 
surement, and the result of (px; + q) is the predicted value of 
the fitted curve, the voltage threshold correction factor fy is 
expressed as follows: 
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This yields an asymmetric coverage interval for fa: 
ETETE, EET According to the GUM, because 
the probability distribution function of fg is unknown, its cov- 
erage interval distribution is assumed to be rectangular by de- 
fault. According to Equation (7), the uncertainty components 
of the voltage threshold correction factors are obtained as fol- 
lows: 
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c. Evaluation of the combined standard uncertainty 
The previous two subsections provide the system dead-time 
and voltage threshold correction factors and their uncertainty 
components, which can be used to obtain the relative uncer- 
tainty of each component. According to the GUM method, 
the relative combined standard uncertainty u( Æ) of the sur- 
face emission rate can be derived from Eq. 1 using the law of 
propagation of uncertainty proposed by Stanga [37]. 


(E + B)? R? RB? 
(E) = RE u?(R) + (E+ By u? fa + T? R?u? (T) + BEIE 7B) (14) 
Here, R and B belong to the Type A evaluation of uncer- 
tainty, both of which were obtained from experimental mea- R 
surements. u(R), u(B), u(fa), and u(r)) are the relative E= =a 1 B (15) 
standard deviations of R, B, fa, and T, respectively. i 2nin2 Rpr+q 


5. Evaluation by the MCM 


The MCM provides a general numerical approach for eval- 
uating the measurement uncertainty for models with any 
number of input quantities that can be characterized by a 
PDF and a single output quantity. In this study, MATLAB 
was used as the Monte Carlo simulation platform to model 
the measurement between the output quantities Y and input 
quantities X,,..., Xy. According to Equations (1), (2), and 
(12), a new measurement function for the surface emission 
rate E can be obtained as follows: 


The probability distributions, best estimates, standard un- 
certainties, and coverage intervals of the propagated input 
quantities X; for each measurement model were analyzed. 
The probability distributions for each input quantity are listed 
in Table 1. 

Finally, the estimated and standard uncertainty of Y were 
calculated from the PDF of Y and the coverage interval of Y 
at a 95% coverage probability (k=2). 


HI. RESULTS 


The order of the experimental measurement was as follows: 
first, the system dead time was measured, then the voltage 


TABLE 1. Information on model propagation inputs for the MCM of measurement. 


Input (X;) Sources of data 


Probability distribution 


R Experimental measurement results Normal 

B Experimental measurement results Normal 

nı Experimental measurement results Normal 

m Experimental measurement results Normal 

ni2 Experimental measurement results Normal 

x Experimental scale readings Rectangular 

p Calculated results Normal 

q Calculated results Normal 
threshold correction factor was determined, followed by the 
count rate and background, and finally, the combined standard ` sao 
uncertainty was calculated. 3000 Curve fitted irom i0-L0Vto WBN) 


1. GUM method 


a. System dead time uncertainty The system dead-time 
measurements were performed using two a-plane sources 
with a count rate of 5000 s71. Ten measurements were taken 
per set of experiments, and the statistical standard deviations 
obtained as the measurement uncertainties are listed in Table 
2: 

The measured values in Table 2 were substituted into Eq. 
(2) to obtain the system dead-time values. According to the 
method described in Section 2.4.1, the half-widths of the in- 
terval of the system dead time, a= 0.2 us, and k = v3, are 
used to obtain the uncertainty of the measurement of the sys- 
tem dead time, as presented in Table 3. 

b. Voltage threshold correction factor uncertainty The 
threshold correction curves of the °’Sr and ”*Tlplane sources 
were measured by increasing the threshold voltage from 
0.040 V. The measured count rate at each set of voltages was 
recorded, and the obtained data are shown in Figures 3 and 4. 
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Fig. 3. The coordinate plots of 8-source (Sr) show eight coordinate 
points, and five coordinate points from 0.10 to 0.18 V were used for 
curve fitting, with the fitted coordinate points and fitted straight line 
marked in red. 
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Fig. 4. The coordinate plots of 8-source ( 7°*T1) show eight coordi- 
nate points, and five coordinate points from 0.10 to 0.18 V were used 
for curve fitting, with the fitted coordinate points and fitted straight 
line marked in red. 


According to the threshold correction curve plots of the 
Sr and ™T] plane sources, the following was found: at 0.04 
to 0.10 V, the curve non-linearity fluctuated greatly, mainly 
because the count rate contained a large number of noise sig- 
nals; at 0.10 to 0.18 V, the correlation coefficient of the fitted 
curves had the largest value of R?, the fitting effect was the 
best, and the normalized values of the count rate (predicted 
value/fitted intercept) were all greater than 90%. The nor- 
malized curve obtained is shown in Figure 5. The above five 
sets of measurements were then entered into Equations (4) 
and (5), and the slopes p and intercepts q of the threshold cor- 
rection curves for the °’Sr and ?“Tiplane sources, as well as 
the standard deviation s(p) of the slopes, the standard devia- 
tion s(q) of the intercepts, the estimated variances s”, and the 
correlation coefficients r(p, q), were obtained and recorded in 
Table 4. 


By substituting the calculation results in Table 4 into Equa- 
tion (11), the uncertainty was evaluated for five sets of data in 
the range of 0.10 to 0.18 V, and the results are listed in Table 
5. 

Both relative standard uncertainties of the two ( sources 
are minimized at the threshold x = 0.10 V. By substituting 
the predicted values at this threshold into Equation (12), the 
voltage threshold correction factors and coverage intervals for 
the two sources were obtained, as shown in Table 6. 


TABLE 2. Dead time measurement results of two-source method. 


Quantity of experiments Measured value (s75 Standard deviation B 


m(Pu) 5229.78 26.10 
no(*'Am) 5403.96 12.61 
nyo(?Pu +”! Am) 10424.83 27.79 
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Fig. 5. Threshold extrapolation straight line for two 8 sources where the counting rates have completed normalization, including error bars 


for each predicted value on the fitted line. 


TABLE 3. System dead time and its uncertainties with the GUM 
method. 


System dead time (s) Standard uncertainty (s) 
3.696 x 107° 0.115 x 10° ° 


The voltage threshold correction factors and uncertainties 
for the two 8 sources were then obtained from Equation (13), 
as listed in Table 7. 

c. Combined standard uncertainty The measured data 
were collated when the a source did not require threshold 
correction by setting its fq value to 1 and its u( fa) value to 0. 
The results are summarized in Table 8. 

By substituting various data from Table 8 into Equation 
(14), the standard and relative standard uncertainties of the 
four plane sources were obtained, as listed in Table 9. 


2. Monte Carlo method (MCM) 


According to the information in Table 1, the input quantity 
X; for the Monte Carlo simulation was defined as follows: 
The input quantities (n1, n2, N12) and their normal distribu- 
tions for the system dead time were obtained from Table 2. 
The standard deviations of p and q were obtained from Equa- 
tions (8) and (9) and were assumed to be normally distributed 
by default. The threshold voltage reading x = 0.10 V had a 
rectangular distribution with a scale component value of 0.02 
V, and its coverage interval was set to [0.08, 0.12]. The mea- 
surements of count rate R and background count rate B in Ta- 
ble 6 were used to obtain the coverage interval of the normal 
distribution of these two quantities using their standard de- 
viations. In the Monte Carlo simulation, the random sample 
size for each input quantity was set to M = 10". The uncer- 


TABLE 4. Parameters of threshold correction curves for two ( sources. 


B sources p q s(p) s(q) 8° r(p,q) 
Sr -989.76 2124.44 61.69 5.09 5.07 -0.98020 
2047] -1160.56 2839.67 52.89 4.36 3.73 -0.98020 


TABLE 5. Standard uncertainties of five sets of predicted values on the threshold correction curves for the two ( sources 


BUG, TOT] 
-Coordinate-Threshold /V 
Standard deviation BT Relative standard deviation Standard deviation (s = Relative standard deviation 

XI 0.10 1.554 0.08% 1.332 0.05% 

x2 0.12 2.618 0.13% 2.244 0.08% 

X3 0.14 3.787 0.19% 3.246 0.12% 

X4 0.16 4.987 0.25% 4.275 0.16% 

X5 0.18 6.200 0.32% 5.315 0.20% 


tainties for the system dead time, {-source voltage threshold 
correction factors, and combined standard uncertainty were 
obtained, and the results are presented in Tables 10, 11, and 
12. 

A normalized histogram of the distribution function of the 
system dead time obtained by the Monte Carlo simulation is 
shown in Figure 6. A normalized histogram of the distribution 
function of the voltage threshold correction factors for the two 
8 sources is shown in Figure 7. A normalized histogram of 
the distribution function of the surface emission rate for the 
four sources is shown in Figure 8. Each histogram contains 
200 bars. 
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Fig. 6. Normalized histogram of the probability distribution function 
(PDF) of the system dead time by the MCM. 


Finally, the value of E, was used to determine whether the 
results of the two methods were consistent within the uncer- 
tainties. The Æ» value is a statistical index used to evalu- 
ate the measurement capability of a laboratory. It is used to 
compare the difference between the laboratory measurement 
results and the reference value, as well as the measurement 
uncertainty. The closer the value of En is to zero, the closer 
the laboratory measurement results are to the reference val- 
ues. If the absolute value of En is less than 1, the laboratory 
measurement results are considered consistent with the refer- 
ence value. If the absolute value of En is greater than or equal 


to 1, it is considered that the laboratory measurement results 
are inconsistent with the reference value. The formula for 
calculating the value of E,, is shown in Equation (16), where 
Ug and Um are the expanded uncertainties (k=2) of the two 
methods. The results are summarized in Table 13. 


(16) 


IV. DISCUSSION 


In the GUM and MCM methods, the system dead-time 
measurement model and voltage threshold correction factor 
measurement methods were the same, but the uncertainty 
evaluation methods were different. Only the 6-plane source 
must complete the voltage threshold correction factor mea- 
surement. The GUM method is based on a Type B evalu- 
ation of the uncertainty process, and the default probability 
distribution of the system dead time and threshold correc- 
tion factors is rectangular. The MCM inputs the measure- 
ments from the system dead time measurement model and the 
voltage threshold correction factor measurement method into 
Equation (15), and the uncertainty results of the system dead 
time and voltage threshold correction factor can be obtained 
through the Monte Carlo procedure. Thus, the final combined 
standard uncertainty can be calculated directly. It should be 
noted that in the MCM, the probability distribution of the 
threshold voltage readings is set to a rectangular distribution 
according to the Type B evaluation of uncertainty; however, 
in reality, the range of the coverage interval of the threshold 
voltage readings is much smaller than the scale component 
value (0.02 V). 

As can be observed from the evaluation results of the GUM 
and MCM methods, the system dead time and voltage thresh- 
old correction factor obtained by the MCM were similar to 
those of the GUM method; however, their measurement un- 
certainties were larger than those of the GUM method. Simi- 
larly, the final measurement estimates of the GUM and MCM 
methods for each source were similar. However, the uncer- 
tainty results of the MCM method were larger than those of 


TABLE 6. Voltage threshold correction factors and coverage intervals at a threshold of 0.10 V. 


B sources q (s71) pxı +q (s71) Standard deviation, px; + q (s71) Coverage interval, pxı + q (s71) Ja 


Coverage interval, fa 


Sr 2124.44 2025.47 1.55 
2047] 2839.67 2723.61 1.33 


(2023.91,2027.02] 
(2722.28,2724.94] 


1.04887 [1.04806,1.04967] 
1.04261 [1.04210,1.04312] 


TABLE 7. Voltage threshold correction factors and their uncertainties for °’Sr and *“Tlof the GUM method 


8 sources fa 


Standard uncertainty 


Sr 1.0489 0.00047 
2047] 1.0426 0.00030 
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Fig. 7. Normalized histograms of the PDF of the voltage threshold correction factors for the two 8 sources (Sr and 2°*T1) by the MCM. 


the GUM method. This can be attributed to two reasons: 
(i) the combined standard uncertainty formula of the GUM 
method ignores the higher-order terms of the Taylor expan- 
sion in the process of partial derivation of the measurement 
model and (ii) the GUM method assumes that the probabil- 
ity distributions of the system dead time and voltage thresh- 
old correction factor have a rectangular distribution, which is 
more complex for both. Both reasons affected the evaluation 
results; however, the author believes that the second reason 
has more weight. 

Compared to the GUM method, the MCM has unique ad- 
vantages: the PDF of the output quantity Y can be output by 
the MCM, and the degree of discreteness of the output value 
can be observed. From the seven PDF histograms, it can be 
directly observed that (i) the probability distribution of the 
system dead time tends to be normally distributed, but the de- 
gree of discreteness is high, (ii) the probability distribution of 
the voltage threshold correction factors of the two 8 sources 
has a trapezoidal distribution characteristic, and (iii) it can be 
intuitively inferred from the PDFs of the four sources that the 
PDFs of ”°Sr have the highest degree of dispersion. 

To improve the measurement method, it is recommended 
to use a digital multichannel analysis measurement system to 
obtain the system dead time directly by measuring the live 
time. Additionally, for a complete measurement system, per- 
formance tests, such as stability and repeatability, should also 


be conducted. 
V. CONCLUSION 


In this study, methods for evaluating the uncertainty of 
a and (6 surface emission rate measurement models based 
on windowless proportional counters were studied using the 
GUM and MCM methods. The measurement methods and 
models of the system dead time and threshold correction fac- 
tor were elaborated in detail, and the uncertainties of the sur- 
face emission rates of the four plane sources were evaluated 
using both the GUM and MCM methods. 


The evaluation results of the GUM method were lower than 
those of the MCM method. However, the MCM method is 
considered to be more objective in its evaluation process and 
can better reflect the distribution characteristics of the surface 
emission rate measurement function through PDF graphs. To 
measure and calibrate the surface emission rate for œ and 8 
particles, the MCM constructs a mathematical model that is 
practical for stable environmental conditions and measure- 
ment methods in laboratories. The GUM method has greater 
flexibility in setting the coverage intervals of the measured 
quantities, whereas the MCM is more suitable for on-site cal- 
ibration. Both methods can be used to validate the evaluation 
results. 


TABLE 8. Values and relative standard uncertainties of the measured surface emission rate of the four plane sources based on the GUM 
method 


Sources E/s-' Bis ' R/s! fa T/s urlB) ur(R) Url fa) Urei(T) 
Tam 5514.70 0.608 5403.96 1 3.696 x 10`" 15.57% 0.23% 0.000% 3.12% 
239Py 5333.47 0.608 5229.78 1 3.696 x 107 15.57% 0.50% 0.000% 3.12% 


Sr 2098.93 39.46 2025.47 1.0489 3.696 x 107 1.15% 0.18% 0.044% 3.12% 
20T] 2827.70 39.46 2723.61 1.0426 3.696 x 107 1.15% 0.32% 0.028% 3.12% 


TABLE 9. Surface emission rate and combined standard uncertainties for the four plane sources of the GUM method 


Sources Surface emission rate (Eg) /s ~! Standard uncertainty (ug)/s ~! Relative standard uncertainty (trei-G) 


Tam 5514.7 13.9 0.25% 
Py 5333.5 27.9 0.52% 
Sr 2098.9 4.5 0.22% 
2047] 2827.7 10.4 0.37% 
35 x105 , , Am-241 oie x108 , , Pu-239 
a b 


Number of Monte Carlo simulations 
Number of Monte Carlo simulations 


-0.02 -0.01 ie} 0.01 0.02 -0.02 -0.01 0 0.01 0.02 
Normalisation of surface emission rate Normalisation of surface emission rate 
x104 Sr-90 
: 


Number of Monte Carlo simulations 
Number of Monte Carlo simulations 


-0.02 -0.01 0 0.01 0.02 -0.02 -0.01 0 0.01 0.02 
Normalisation of surface emission rate Normalisation of surface emission rate 


Fig. 8. Normalized histogram of the PDF of the surface emission rate for the four sources (Am, **°Pu, Sr, T1) by the MCM . 


TABLE 10. System dead time and its uncertainty with the MCM 


Dead time (s) Standard uncertainty (s) 
3.694 x 1078 0.697 x 107° 


TABLE 11. Voltage threshold correction factors and their uncertain- 
ties for °°Sr and *'Tlof the MCM 

B sources fa Standard uncertainty 

Sr 1.0489 0.0047 

MT] 1.0426 0.0055 
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TABLE 12. Surface emission rate and combined standard uncertainties for the four plane sources of the MCM 


Sources Surface emission rate (Em) [s7] Standard uncertainty (u Myst Relative standard uncertainty (Urei-c) 


Tam 5513.5 24.9 0.45% 

Py 5332.3 33.6 0.63% 

Sr 2101.0 14.6 0.69% 

24T] 2829.1 18.6 0.66% 
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